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ABSTRACT 


The behavior of shallow spherical caps subjected to 
axisymmetric and nearly axisymmetric step pressure loads is 
examined using a digital computer program for the geo- 
metrically nonlinear analysis of arbitrarily loaded shells 
of revolution. A criterion for dynamic buckling under the 
nearly axisymmetric load is proposed and critical buckling 
pressures are determined for both loading conditions for a 
large range of cap sizes. The axisymmetric results either 
compare well with, or are somewhat higher than, previously 
published theoreticai results. For most of the cap sizes 
the buckiing pressures Irom the asymmetric analysis are 
new results. A numerical convergence study is conducted 


to verify the accuracy of the results. 





TABLE OF CONTENTS 


i. INTRODUCTION ------------------------------------ 6 
II. DESCRIPTION OF THE COMPUTER PROGRAM ------------- 8 
III. PROBLEM DESCRIPTION ---------------------------- - ll 
IV. BUCKLING CRITERIA ------------------------------- ee 

A. AXISYMMETRIC ANALYSIS ---------------------- - 13 

B. ASYMMETRIC ANALYSIS ------------------------- 14 
,. PV ORIN Sg ee ere 16 

A. AXISYMMETRIC ANALYSIS ----------------------- 16 

B. ASYMMETRIC ANALYSIS ------------------------- 17 
VI. DISCUSSION -------------------------------------- 19 
Vil. WNUMERICAL CONVDORGENCE sTuDY --------------------- my 
VIII. SUMMARY AND CONCLUDIONS ------------------------- 24 
FIGURES AND TABLES ------------------------------------ 26 
LIST OF REFERENCES ------------------------------------ 34 
INITIAL DISTRIBUTION LIST ----------------------------- 35 
FORM DD 1473 ------------------------------------------ 36 





K (n) 


U, V, W 


TABLE, OF SYMBOLS 


the modulus of elasticity of the shell 


the rise of the spherical cap at the pole 
the thickness of the shell 


the square matrix associated with the linear portion 
of the governing equations 


the mass density of the shell 

the meridional bending moment per unit length 

the Fourier index 

a nondimensional applied load 

the nondimensional critical pressure 

the classical buckling pressure of a complete sphere 


a column 


matrix containing the coefficients of the 
nth term 


in the series expansion of the applied load 


a column 
temms ain 


matrix containing the inertial and nonlinear 
the nth set of equations 


the normal distance from the axis of revolution to 
the surface of the cap 


the normal distance from the axis to the cap in the 
base plane; the maximum value of r 


PaCmaacwInOL CUmvarlre in the s ande@6 directions 
respectively 


the meridional distance along the surface of the 
shell 


the nondimensional time 
the time 
a Gweterence time 


the displacements in the s, 9 and ¢ directions 


respectively 





(n) 


aenongumensrvonal measure of the volume of the shell 
deformation 


the peak in the time history of the parameter V 


the displacement in the Zt direction in the nth 
harmonic 


the nondimensional time increment 


the nondimensional parameter governing the magnitude 
of the load applied in the asymmetric harmonics 


the CoOerainace Normal to the surface of the shell 


the circumferential angle measured about the axis 
of revolution 


a nondimensional poouSeene parameter used to describe 
the spherical cap 


Poisson's ratio 


the normal distance from the base plane to the middle 
surface of the undeformed cap 





i. PENTRODUCTION 


The critical pressures for the static axisymmetric snap 
buckling of uniformly loaded clamped shallow spherical caps 
were presented by Huang [1] in 1964. Results for the 
critical pressures for dynamic axisymmetric buckling of 
clamped shallow spherical caps subjected to a step pressure 
load have been recently given by Huang [2], by Stephens 
Saas Fulton [3], by Stricklin et al. [4] and by Ball [5]. 
The results for the axisymmetric dynamic buckling indicate 
that the critical dynamic pressure varies from 40% to 75% 
@= the Critical pressure for static,axisymmetzrac buckling, 
@epeaaing upon the size of the cepe i=in nis Static analy sus, 
Huang [1] shows that asymmetric buckling modes exist at 
pressures below the critical pressure for static axisymmetric 
Snap buckling. These asymmetric modes are unstable and 
cause the cap to buckle at the minimum bifurcation buckling 
pressure. If imperfections are present the buckling will 
occur at pressures below the minimum bifurcation buckling 
pressure. 

This thesis considers the possibility that the same 
Situation may exist for the dynamically loaded cap, i1.e., 
the presence of asymmetric motion may reduce the critical 
dynamic buckling pressure. A digital computer program for 
the geometrically nonlinear analysis of arbitrarily loaded 


shells of revolution [5] is used to determine the response 





of several caps when subjected to axisymmetric and nearly 
axisymmetric step pressure loads. It was anticipated that 
the presence of the very small asymmetric portion of the 
load would initiate asymmetric motion and that this motion 
might grow large and significantly influence the axi- 
symmetric response. If the axisymmetric response is not 
influenced by the asymmetric motion, the buckling is 
referred to as direct snapping; if the asymmetric motion 
reduces the critical pressure, the buckling is called 
indirect snapping [6]. Stricklin et al. [4] have also con- 
sidered the dynamic asymmetric buckling problem and have 
presented a range of critical pressures for three cap sizes. 
The objective of the thesis is to determine more precisely 


A “ 5 = = a a er " gs Oe Se ee = st ~~ : ~ 
Brceeiertacei Gynammc pressure for 2 marge Tange @i sep Sazes. 





me. WeSoCRte wae OF THE COMPUTER PROGRAM 


The computer program used in this study is capable of 
analyzing thin, linearly elastic shells of revolution sub- 
jected to arbitrary load and temperature distributions. 
The governing partial differential equations are based upon 
Sanders’ nonlinear thin shell theory for the condition of 
small strains and moderately small rotations [7]. The 
inplane and normal inertial forces are accounted for, but 
the rotary inertial terms are neglected. The set of 
governing nonlinear partial differential equations is 
reduced to an infinite number of sets of four second-order 
Gdiffererntiai equaticws in the meridional end time coardi~ 
nates by expanding all dependent variables in a sine or 
cosine series in terms of the circumferential coordinate 6@. 
The sets are uncoupled by using appropriate trigonometric 
identities and by treating the nonlinear coupling terms as 
pseudo loads. The meridional derivatives are replaced by 
the conventional central finite difference approximations, 
and the displacement accelerations are approximated by the 
implicit Houbolt backward differencing scheme [8]. 

This leads to sets of algebraic equations in terms of 
the dependent variables and the Fourier index, n, of the 


form 


er re gy eta 4 (1) 





The square Imatrix (x, ()] is a banded matrix associated with 
the linear portion of the governing equations. The elements 
in the column matrix (z'™) are the coefficients of the nth 
term in the series expansions of U, V, W, and Me at each 
finite difference station where U, V. and W are the displace- 
ments in the meridional, circumferential and normal direc- 
tions respectively, and M. is the meridional bending moment 
per unit length. The applied loads in the ace harmonic are 


a 


contained in the colum matrix {gq and the inertial and 
nonlinear terms are contained in the column matrix igen. 
At each time step, an estimate of the solution is 
obtained by extrapolation from the solutions at the previous 
time steps. An elimination procedure is used to repeatedly 
sOlve each set of aluebraic equations until the maximum 


difference between two consecutive solutions {2Z is 


(n) 5 
ii 2,..., 1S Less than the convergence criterion, which 

is a specified percentage of the maximum value of U, V, W, 
and M. in the nt? parmonic. A time history of the shell 
response is obtained by incrementing the time in equal steps 
until either a prescribed maximum number of time steps have 
been taken or the solution fails to converge after a speci- 
fied number of iterations. The time interval is usually 

so small only a few iterations are required since the 
difference between the estimated solution and computed 
solution is generally small. However, when the shell 
becomes dynamically unstable, the solution may not converge, 


even with many iterations. 





The program is dimensioned such that the product of the 
number of meridional stations and the number of harmonics 
must be less than 201, and the maximum number of harmonics 
must be less than ll. The program occupies approximately 
155,000 bytes on the IBM 360/67 Digital Computer, anda 
typical execution time using 40 stations and 2 harmonics 
with 3000 time steps was 12 minutes using the FORTRAN IV, 
LEVEL H Compiler, with OPT = 2. A copy of the program for 
a static analysis and USERS MANUAL can be obtained through 
COSMIC (M70-10098,LAR-10736). The program with the dynamic 
analysis capability will be available in the latter half of 


ie. 
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ITis SeseeGeEM DESEGRIPTION 


The geometry of the shallow spherical cap can be speci- 


fied by the nondimensional parameter i, where 
1 1. 
dh = 2[3(1-v*) ]*(H/h)* (2) 


and H is the rise of the cap, h is the thickness and v is 
Poisson's ratio. The mass density of the cap is m. For 


ee shells considered here: 


Radii of curvature, Ros Ro = 250 in. 
Thickness, h => Ore 2 5 dei 
Modules of Elasticity, E = Mee s10° psi 
Moascan's tateic, v POs 


The classical buckling pressure of a complete sphere do is 


given by 
qo = OE tnr,)* (3) 
ey) 


The nearly axisymmetric step pressure load is given by 


{q')} = eels 
(4) 
iq} = Pq.e {1}, i Sp ae 
(n) 


in which P is a nondimensional load parameter, € is a 


peel number and T is time. The shell is stress free and 


mMetioniess prior to the application of the load. For the 


axisymmetric analysis, only the n = 0 harmonic is considered. 


lek 





In the asymmetric analysis various harmonics were excited 


(n) 


either singly or in groups through the parameter € =0.0001. 
The nondimensional time is denoted by t, where 
aanayn (5) 
and 
T = (R_?m/E)® (6) 
O S 


For this study, m was selected such that Te was equal to 
one. For most of the runs, the non-dimensional time incre- 
ment dt was taken as 0.05 and the maximum response time was 
150. The necessity for the long response time was noted 


mieRererencer 3. 
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iV BUGKIOING ‘CRETERIA 


Criteria for dynamic buckling of shells are not as well 
defined as they have been for the case of static buckling. 
Quoting from Stoker [9], ". . . stability in dynamics in an 
idea that means different things to different people, .. ." 


and "Criteria for stability that are logically sound, and 


also practically applicable, are not easy to come by either." 


A. AXISYMMETRIC ANALYSIS 

The criterion most often used for dynamic axisymmetric 
buckling of the spherical cap is based on the variation of 
Veen the peak in the time history of a nondimensional 
average displacement parameter V, with respect to P, the 
amplitude of the applied load [2, 3]. The parameter visa 
measure of the volume of the shell deformation and has been 


defined in Reference 3 as 


vo 
f rw (9) 6le 
v= —_—_—_ CT) 
Oo 
juan’ Sel Oa 


(> 


in which r is the normal distance from the axis to the cap, 
ry is the maximum value of r, and & is the vertical distance 
from the base plane to the undeformed cap middle surface. 


For loads below the dynamic buckling load, where the non- 


linearities are small, the relationship between V and P is 


les 





approximately linear. There may exist, however, a certain 
value of P at which a very small increase in P produces a 
very large increase in SU net This load is defined in 


References 2 and 3 as the buckling load. 


B. ASYMMETRIC ANALYSIS 

At the present time no criterion for the dynamic buckling 
of shallow spherical caps under nearly axisymmetric loading 
conditions has been proposed. The results of this study 
suggest that the criterion for buckling in the asymmetric 
a@alvysis #s identical to the criterion for buckling In gene 
axisymmetric analysis, 1.ée., buckling is defined to occur 
when a small change in the nearly axisymmetric applied load 
causes a large change in Veo? where V is given by equation 
(7). This criterion is based upon the argument that the 
axisymmetric deformation is the only true measure Oi 
buckling. The cap can suffer large displacements el ele 
asymmetric harmonics, but these displacements in themselves 
do not represent a catastrophic buckling failure because 
the transverse displacements at the pole and base of the 
cap are zero in the asymmetric harmonics. Thus, the shell 
will have snapped through, or buckled, only when V, the 
volume of deformation, becomes large. The onset Of buck ling 
may occur at loads smaller than those that cause buckling 
in the axisymmetric analysis because of the influence of 
motion in the asymmetric harmonics on the axisymmetric 


response. 
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For convenience, the response in each asymmetric harmonic 
will be measured using equation (7) with w (0) replaced with 
wir) The parameter V for the asymmetric harmonics does not 
represent a volume of deformation as it does for the axi- 


symmetric harmonic. It can, however, bestised@ “covindacate 


the relative excitation of the asymmetric harmonics. 


5 





VV. RESULTS 


Po hoe So YMMETRIC ANALYSIS 

Plots of V versus t for two values of P are given in 
Figures 1 through 3 for X} = 5, 6 and 9 respectively. In 
Figures 1 and 3 one value of P is slightly above the criti- 
cal pressure. In Figure 2 both values of P are slightly 
above the critical pressure. Note that the shells \ = 5 
and 6 buckle considerably earlier than the shell A = 9. 


Mae shells X} = 8 and 10 also buckled relatively early, while 


the shells } = 7, 7.5, ll, 12 and 13 responded in a manner 
Similar to the shell X = 9 and oscillated many cycles before 
buckling. When the applied prassure was somewhat Larger 


than the critical pressure, the time to buckle for most of 
the shells decreased, as illustrated in Figure 2. However, 
the shell A = 13 took longer to buckle for a pressure 


approximately 10% larger than the critical pressure. 


Replet of VMAX 


eee. ., L3 and 7.5.2 Note the distumet break in Ehe rela 


tionship between P and Ve for all X except A = a, In 


versus P is given in Figure 4 for ir = 4, 


every case the point after the break was accompanied by a 


Diiie results for the shells A = 5, 8, and IY were first 
presented in Reference 5. 

2 The Sem A = 4 did not exhibit a distinct break in 
the P - Vyay relationship. Hence, there is no well defined 
critical pressure. Nevertheless, there is a large change 
in Vyay over a relatively small range of pressure. A simi- 
Har Satwatien exists for the static buckling of shells with 
small i. 
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convergence failure. Hence, larger values of Mae would 
have occurred if the program could have continued. The 
value of P at the break is defined as the critical pressure, 
Porir’ and is given in Table I and in Figure 5 for the 
values of A considered. Figure 5 also shows the results 
presented by Huang [2], by Stephens and Fulton [3] and by 
Stricklin et al. [4]. The results for static axisymmetric 
buckling presented by Huang [1] and by Stilwell [10], who 


used the same program as the present study, are also shown 


in Figure 5 for purposes of comparison. 


B. ASYMMETRIC ANALYSIS 

Due to the significant increase in computer time required 
to examine several harmonics simultaneously, only a few of 
the harmonics in the vicinity of the critical harmonic for 
static buckling were excited. To verify the validity of 
this procedure three shells were examined using additional 
harmonics. These results will be presented in another 
Section. | 

As an example of a typical response curve, a plot of 
both wh) ry mide at © = 0.513 Ly for the harmonics n = 0 
and 2 is given in Figure 6 for the shell A = 6 with 
Me=sOu@ieo. The critical harmonic for static buckling of 
Cieemeamell 15 n = 2. Note in Figure 6 the agreement 


(1) 3 and V in both harmonics. Also shown in 


between W 
Figure 6 is V for n = 0 and P = 0.566 without the presence 


of the asymmetric harmonics. 
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Similar runs were made for =o Orme, 4 6 ous ee omana «7.5 
with several values of P. Only the shells A = 6, 7, 7.5, 
ll and 12 were influenced by the asymmetric motion. The 
other shells were not influenced by the motion in the 
asymmetric harmonics and buckled at the critical pressure 
from the axisymmetric analysis. The plots of Vinee Teor 
n = 0 versus P are shown in Figure 7 for the five shells. 
Note the distinct break in the Viapeats - P relationship. 
Furthermore, the point after the break was accompanied by 
a convergence failure. However, the point at P = .52 for 
the shell X} = 6 was a converged solution. 

The critical pressures for all the shells examined in 
the asymmetric analysis are given in Table II and in Figure 
metre range of Cfitical pressures Peporvea by BStrieklin 
et al. [4] are shown in Figure 8 by bars encompassing the 
span of values for which they found appreciable excitation 
of asymmetrical displacement coefficients. The experimental 


results for two shells with 4A = 7.5 presented by Lock 


et al. [6] are also given in Figure 8. 


lee 





Wie DISCUSSION 


The major results of this study are the critical pres- 
sures from the asymmetric analysis given in Figure 8. These 
results show that when the critical pressure from the axi- 
symmetric analysis is larger than approximately 0.52, the 
Shell will buckle at P = 0.52 due to the presence of 
mee2on in the critical asymmetric harmonic. In this 
analysis, the motion in each asymmetric harmonic was caused 
by a load of .01% of the step pressure load. An analysis 
that uses very small asymmetric imperfections in the shape 
of the cap as the triggering mechanism should give the same 
FesuiLts. 

The two experimental buckling pressures presented in 
Reference 6 for the shell A = 7.5 are below the value pre- 
dicted in this study. The discrepancy is probably due to 
the imperfection sensitivity of the spherical cap. Accord- 
ing to Reference 6, the response of the two shell models 
was Significantly different. One shell buckled shortly 
after the load was applied; the other shell vibrated for 
about nine cycles prior to snapping. The reason for this 
can be seen in Figure 8. The shell A = 7.5 is almost a 
borderline case and the buckling may or may not be influenced 
by asymmetric motion, depending upon the exact properties 
of the shell. In Reference 6, direct snapping was 


characterized by a rapid growth in the axisymmetric response 
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with little oscillation. An example of this type of 
behavior is shown in Figure 1. On the other hand, Figure 

3 with P = 0.44 illustrates a different kind of direct 
snapping behavior since many oscillations occur prior to 
buckling. However, an additional increase in the load 
caused most of the shells to buckle much sooner as illus- 
trated in Figure 2. The a snapping discussed in 
Reference 6 was characterized by several oscillations in 
the axisymmetric harmonic with a steady growth in an asym- 
metric harmonic until snapping occurred in the axisymmetric 
mode. The response shown in Figure 6 is an illustration 
of this type of behavior. Thus, a shell may oscillate 
several cycles before buckling if either the load is equal 
Manos Siightly above the criticel Loa@ of the si@éf8l Wuckles 
due to asymmetric motion. Hence, the shell A = 7.5 could 


exhibit both direct snapping and indirect snapping. 
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VII. NUMERICAL CONVERGENCE STUDY 


The numerical parameters used in the analysis could, if 
not selected properly, lead to inaccurate solutions. Conse- 
quently, the influence of the time step, the convergence 
Griterion, the number of finite difference stations along 
the meridian, double precision versus single precision, and 
the assumed critical dynamic asymmetric harmonic was 
examined for several values of i. 

The time step was reduced from 0.05 to 0.025 for the 
axisymmetric analysis of the shells } = 5 and ll. In 
References 2 and 3, 6t = 0.02 and 0.10 respectively. For 
the shell A = 5, the critical load was 0.46, aS before; 
but the solution at P = 0.49, which displayed a sharp rise 
in Vuaxs 


smaller time step allowed the solution to enter the post- 


did not experience a convergence failure. The 


buckling region. For the shell A = ll, Porit = 0.575, and 
the solution failed to converge at P = 0.58, as before. 
The results presented in Figures 1 - 8 were obtained 


using a convergence criterion of .01. To see if this was 
adequate, the critical axisymmetric pressure of the shell 

AX = 10 was determined using several smaller values. 
Decreasing the value of the criterion to 0.001 and to 

0.0001 resulted only in requiring more iterations (a maximum 
of two were required at 0.01, three were required at 0.001 


and four were required at 0.0001) with no change in the value 


Px 





Or VMAX and PoRItT’ The computer execution time increased 
from approximately six minutes for 0.01 to approximately 
twelve minutes for 0.001 and eighteen minutes for 0.0001. 
The number of difference stations used to model the 
shell was forty. This agreed with number found to be 
wieeacient for the static analysis for A < 10 [10]. Twenty- 
Six stations were used in Reference 3 and Reference 2 repor- 
ted using 4A + 1 stations. The number of stations was 
increased to eighty for the axisymmetric analysis of the 
shell } = 11. No significant change was found in the 
results. 
The computations for all the shells were made using 
the single precision mode of the IBM 360/67. Results for 
the shells A = 6 ana iO were aiss cktained using the double 
precision mode. The double precision arithmetic produced 
no appreciable change in the critical axisymmetric pressures. 
In order to establish the validity of the assumption 
that the critical asymmetric harmonic is in the vicinity 
of the critical harmonic for static buckling, the solution 
for the shell A = 6 was obtained using the four harmonics 
me—- OO, Ll, 2 and 3. The critical harmonic for static buck- 
ling is n = 2. The results showed that the response in the 
harmonics n = 1 and 3 was relatively insignificant, but, 
as shown in Figure 6, V for the harmonic n = 2 was of the 
same order of magnitude as V for the axisymmetric harmonic. 
Hence, the critical asymmetric static harmonic appears to 


also be the critical asymmetric dynamic harmonic. A check 
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on the critical asymmetric harmonic was also made for the 
[am XA = b2, The critical static harmonic is n= 7. Com- 
Pimeedons Of harmonics from n = 4 to n = 9 were excited. 


Although there was no appreciable change in the critical 


pressure, harmonics n = 5 and 6, not 7, showed the largest 
motion. The shell A = 9 with the critical asymmetric static 
harmonics n = 4 and 5, was examined using combinations of 
harmonic n= 1ton= 5. None of these harmonics was 


appreciable in comparison with the axisymmetric motion, nor 


did they have any effect on the critical pressure. 








VIII. SUMMARY AND CONCLUSIONS 


A digital computer program for the geometrically nonlinear 
analysis of arbitrarily loaded shells of revolution has been 
used to determine the minimum buckling load of clamped shal- 
low spherical caps subjected to axisymmetric and nearly 
axisymmetric step pressure loads. For both loading condi- 
tions, the critical pressure was defined to be that pressure 
which when given a small increase caused a large change in 
the volume of deformation. In the axisymmetric analysis, 
the critical pressure was determined for the shells i = 5, 
6,.-.-, 13 and 7.5. In general, the results either agree 
well or are somewhat nigher than previousiy pubiishec 
results. In the asymmetric analysis, motion was introduced 
in several of the harmonics in the vicinity of the critical 
harmonic for static buckling. Those shells with a critical 
pressure from the axisymmetric analysis greater than about 
52% of the classical buckling pressure of a complete sphere 
weeewoe 7, 7.5, ll and 12) buckled at a pressure ratio of 
approximately 0.52 due to the asymmetric motion. The other 
shells (A = 5, 8, 9, 10 and 13) buckled at pressure ratios 
below 0.52 and were not affected by the asymmetric motion. 

These results were obtained using a perfect shell of 
revolution with a very small load in each asymmetric har- 


monic. The critical pressures of actual imperfect shells 
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will be less than the pressures presented here. How much 
less depends upon the nature of the initial imperfections 


and is a subject for future study. 
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